NASA NASA -TP-2351 19840020681 

Technical 

Paper 

2351 

August 1984 

A Numerical Simulation 
of Three-Dimensional 
Flow in an Adaptive 
Wall Wind Tunnel 

J. P. Mendoza 







NASA 

Technical 

Paper 

2351 

1984 


A Numerical Simulation 
of Three-Dimensional 
Flow in an Adaptive 
Wall Wind Tunnel 


J. P. Mendoza 

Ames Research Center 
Moffett Field , California 


NASA 

National Aeronautics 
and Space Administration 


Scientific and Technical 
Information Branch 




SYMBOLS 


a‘>* 

$ 

(a«T l 


b 


C L 




MSE 


w 


influence coefficient, A w l /Aw j p influence 
coefficient matrix 

element of the inverse matrix 

inverse of matrix a l 1 


wing span, cm 
lift coefficient 

mean aerodynamic chord, cm 
free-stream Mach number 


mean square error, 



vertical velocity normalized by the free- 
stream velocity 


Wp normalized vertical velocity through panel 

centroid 


x,y,z 

a 

Aw 
A w* 


A J p 

<P 


Cartesian coordinates 

angle of attack, deg 

difference in vertical velocities 

difference in vertical velocities at the ith field 
level control point 

difference in vertical velocities at the jth panel 
centroid 

perturbation potential function 


4> x ,(l)y,<j> z d<t>lbx, btplby, d<t>/dz 

<t> xx dyy^zz d 2 < t>l dx2 > d 2 ( Pl d y 2 ’ <> 2 <t>ldz 2 




SUMMARY 


Numerical simulations of three-dimensional flows in a prototype adaptive wall wind 
tunnel were conducted at the Mach number of 0.6 to investigate : (1) wind-tunnel wall 
interference, (2) active streamline control by varying air removal or injection along the 
walls , and (3) to develop a method for establishing wall boundary conditions for 
interference-free flows . It was found that wind-tunnel wall interference could be con- 
trolled by using only the vertical velocity components. For the configuration tested , it was 
found that interference-free flow with solid sidewalls can be approximated by using only 
floor and ceiling blowing/ suction. 


INTRODUCTION 


The adverse effects of wind-tunnel wall interference have 
always been considered an important factor in wind-tunnel 
testing. An early scheme to reduce wall interference in high- 
speed flow (refs. 1 and 2) involved the use of inserts to con- 
form to a theoretically determined interference-free wall 
geometry. Later (ref. 3), jacks at several stations along the 
floor and ceiling were used to shape flexible walls to theoreti- 
cally determined contours. These early efforts to reduce 
wall interference were soon discontinued because the calcu- 
lations that were required to obtain the required wall con- 
tours were too laborious. Experimental studies in wind- 
tunnel wall interference have recently resumed. Present day 
investigations differ from those of the early forties princi- 
pally in the use of the adaptive-wall method introduced by 
Ferri and Baronti (ref. 4) and Sears (ref. 5). The concept, 
which was made practicable by the use of high-speed com- 
puters, is to measure two different flow quantities (such as 
the axial and vertical components of the velocity along a 
given line or plane) in the wind tunnel. One of the flow 
quantities is used to compute the interference-free values of 
the other quantity by using certain functional relationships 
that exist between the two quantities. For example, in one 
application (ref. 5) measured axial velocities were used to 
compute vertical velocities, the differences between the com- 
puted and measured vertical velocities being a measure of the 
wind-tunnel wall interference. Wall boundary conditions are 
adjusted and wall interference is re-assessed. This process is 
repeated until the differences between the measured and 
computed velocities are sufficiently small so that the flow is 
essentially interference free. 

There are two main methods currently being used to 
adjust the wall boundary conditions. The streamlined -wall 
method (refs. 6 to 1 1) uses a series of jacks for control of the 
wall contour. The second method attempts to effect 


streamline changes within the flow by actual blowing or suc- 
tion at the ventilated walls of the tunnel (refs. 12 to 19). The 
adaptive-wall studies in the Aerodynamics Division of Ames 
Research Center have concentrated on the latter approach 
(refs. 20 and 21). 

The method adopted to assess wall interference (ref. 22) 
uses one measured flow quantity at two surfaces to assess 
residual wall interference. The method has been successfully 
demonstrated in two dimensions and an extension to three- 
dimensional flow has been reported in reference 21 . 

Because of the many parameters that are involved in the 
design of three-dimensional adaptive-wall wind tunnels, it is 
important to simulate the iterative process that would be 
followed in the wind tunnel as faithfully as is possible. In 
reference 23, two-dimensional and axisymmetric computer 
simulations of the adaptive-wall wind tunnel were used to 
demonstrate a “one-step iteration” to interference-free flow. 
In this paper, a new computer simulation is described and is 
used to drive the wind-tunnel flow to an interference-free 
flow condition without any a priori knowledge of the wall 
boundary conditions. 

A theoretical investigation is conducted here in numerical 
flow simulations on a semispan wing that is mounted on the 
sidewall of an adaptive-wall wind tunnel. The flow simula- 
tions were conducted at = 0.6 at two angles of attack, 
1° and 2°. The computer code (outer flow solver) that was 
developed for the present study will be used to develop a 
numerical adaptive wall simulator for configuration design 
studies of a new adaptive-wall wind tunnel. 

METHOD 


Implementation of the adaptive-wall wind-tunnel concept 
(refs. 4 and 5) requires: (1) measurement of flow variables; 
(2) assessment of wind-tunnel wall interference; and 



(3) adjustment of wall boundary conditions. These steps are 
repeated until interference-free flow, as computed in step (2) 
above, is indicated. The operation of an adaptive-wall wind 
tunnel can be represented by the flow diagram shown in fig- 
ure 1 . Measurements of the flow variables are made in the 
wind tunnel at two imaginary planar surfaces surrounding the 
model. The measurements on the planar surface that is 
closest to the model (the source level) are used as inputs to 
the wall-interference assessment program from which the 
values for Aw are computed, the values being the differences 
between the vertical velocities that are computed by the 
outer flow solver and those at the planar surface farthest 
from the model (the field level). If the values for Aw satisfy 
the convergence criteria, the operation is terminated; other- 
wise, they are used as inputs to the influence coefficient 
program to determine a new set of wall boundary conditions 
(Aw^) for the wind tunnel. The blocks labeled “wind-tunnel 
flow field,” “wall-interference assessment,” and “influence 
coefficient” will be discussed in more detail. 


Wind-Tunnel Flow Field 

The wind-tunnel flow field is modeled by a particular 
solution to Laplace’s equation 

$xx + ^yy + ^ = 0 0 ) 

Although a number of computer programs are available to 
solve the potential flow equation with wind-tunnel boundary 
conditions, the Douglas-Neumann Potential Flow Computer 
Program (ref. 24) was selected for this study because it has 
the following features: (1) the program can be applied to 
complex configurations; (2) the user can specify velocities 
on the surface of the configuration being studied; and (3) the 
user can specify points in the flow field where the flow quan- 
tities are to be computed. The computer program is a “panel 
code” that uses the Gothert transformation to account for 
the effects of compressibility (Mach number). 

The wind tunnel paneling scheme is shown in figure 2. 
The wind tunnel is 53-cm long, 25-cm wide, and 13-cm high. 
These dimensions model the Ames adaptive-wall wind tunnel 
used in the experiments reported in references 20 and 21. 
The surfaces of the wind tunnel were represented by 
168 panels. The floor and ceiling surfaces were each repre- 
sented by 36 panels (three panels in the lateral direction and 
twelve panels in the longitudinal direction) and each of the 
sidewalls were represented by 48 panels (four panels in the 
vertical direction and twelve panels in the longitudinal direc- 
tion). The wing (fig. 3) was supported from one sidewall 
(which serves as a reflecting plane) with its leading edge at 
two chord lengths from the upstream-end of the test section. 
It was represented by 200 panels on each of the upper and 
lower surfaces. The coordinate system used in this investiga- 
tion is shown in figure 4. 


Although three velocity components may be specified at 
any of the panel centroids, only velocities normal to the 
floor and ceiling panels were introduced as boundary condi- 
tions in the current study. (These velocities were meant to 
simulate actual blowing or suction at the floor and ceiling of 
the wind tunnel.) It was assumed that solid sidewalls were 
used in this investigation. Figure 5 shows a typical vertical 
velocity distribution (at 2 y/b = 0/ 71, z/c = 0.63) induced by 
blowing through a singles point in the ceiling (tunnel empty) 
at M oo = 0.6 where the blowing was simulated by a source at 
the panel centroid. 


Wall Interference Assessment 

The method of assessing wind-tunnel wall interference in 
this study follows the one described in reference 22 which 
used only the vertical velocity component. In solving the 
‘ ‘wing-in -the -wind -tunnel” problem, vertical velocities were 
computed at control points on two different imaginary 
planar surfaces (the source and field levels) that enclosed the 
model (fig. 6). The field level control points are at 
x/c =-0.69, -0.46,-0.23,0., 0.23,0.46, 0.69, 0.92, and 1.15. 
The point x/c = 0 is at cl 4. The lateral spacing of the control 
points is 2 y/b = 0.24, 0.71, 1.18, and 1.30. The vertical 
spacing is z/c = ±0.52 and 0.0. Using the vertical velocities on 
the surface closest to the model (source level) the vertical 
velocities on the other surface (the field level) were calcu- 
lated by solving a second potential flow problem (the “outer- 
flow solution”). Specifically, the vertical velocities, (j> z = w, 
on the source level together with the far-field boundary 
conditions, <j> x = <p y = <p z = 0, are used to numerically solve 
the compressible small-disturbance equation 

(1 - xx + (j)yy + <p zz - 0 ( 2 ) 

The solution to this second potential flow problem (the 
outer flow solution) is used to compute the vertical velocities 
on the field level surface. The differences between the pre- 
viously computed vertical velocities on the field level surfaces 
and those computed using the outer flow solver are a mea- 
sure of wall interference. When these differences are suffi- 
ciently small, the flow can be considered interference free. In 
this study the mean square error (MSE) based on the differ- 
ences between the vertical velocities that are computed by 
the method in reference 24 and those that are computed by 
the outer flow solver was used as a measure of the residual 
interference. When this parameter is less than or equal to a 
given minimum value, then the flow is considered to be inter- 
ference free. The minimum values (see appendix) are 
0.493X 10 -3 for the upper and lower field level surfaces 
{z/c = 0.52) and 1 .099X 10 -3 for the station that is outboard 
of the wing tip (z/c - 0). (In this study all the velocities have 
been normalized by the free-stream velocity.) The minimum 
values for the MSEs are discussed further in the appendix; 
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Influence Coefficients 

The method of influence coefficients used in the present 
study follows the method first used experimentally in refer- 
ence 20. Basically, an efficient method is needed to deter- 
mine the required increment in the boundary conditions at 
each of the floor and ceiling panels when the increment in 
vertical velocities is given at each of the control points on the 
field level surface. Specifically, the change in the velocity, 
Aw*, at the zth field level control point is equal to a linear 
combination of the velocities, A w|>, through the floor and 
ceiling panels. This is expressed by the equation 
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Here, the a lJ are the influence coefficients which are assumed 
invariant with each iteration. The procedure for determining 
the a l f is as follows: velocities (A w|>) were specified at one 
pair (in this case the jth pair) of opposite -facing floor and 
ceiling panels. Flow entered the test section (blowing) at the 
floor and exited (suction) through the ceiling at a single pair 
(the yth pair) of floor and ceiling panels. The velocities, 
Aw 1 , Aw 2 , . . . , Aw^ that were induced at the field level 
control points were used to compute influence coefficients, 
a 1 ! = AwVAw^, column by column. Assuming that there is 
an equal number of control points and panels, the a 1 ^ -matrix 
in equation (3) is square and may be readily inverted. Assum- 
ing that an inverse matrix {aV)~ l exists, then the following 
equation can be obtained by multiplying both sides of 
equation (3) by (fl#)" 1 . 


1 


~b 11 

l b 12 

A w 2 

= 

b 21 

b 22 

Aw™ 


b mi 

b m2 


b im 
b 2m 

Mm 



Aw 1 


Aw 2 

X 

Aw m 


(4) 


The left-hand side of this equation represents^the vectors of 
the required wall panel velocities and b l 1 = ( 1 .When the 
values for a lJ were computed, it was found that the strongest 
influences were produced by the blowing or suction on the 
nearest wall panel. For example, the field level control points 
that were closest to the ceding were quite insensitive to the 
flows through the floor panels. This effectively decoupled 
the control of the upper and lower halves of the test section. 


A similar result was reported in the experimental investiga- 
tion in reference 20. It was also found that the subset of 
panels shown in figure 7 produced the strongest influence on 
the velocities at the field level control points depicted by the 
crosses. The induced velocities from the other panels were so 
small that they were neglected. The final influence coeffi- 
cient matrix was an 18 by 18 array. 


RESULTS AND DISCUSSION 


The first part of the present investigation was to establish 
procedures for effecting interference-free flows in the wind 
tunnel. This was done with the wing at 2° angle of attack. 
The second part was to apply those procedures to the same 
wing at 1° angle of attack. The Douglas-Neumann code failed 
when applied to the a = 3° case. 

A natural starting point for the beginning of the adaptive- 
wall wind-tunnel operation would be the computation of the 
solid-wall wind-tunnel flow characteristics. Figure 8 shows 
the vertical velocity distributions that are computed at the 
field level control points and those that are computed by 
using the outer flow solver. Note the large wall interference 
effects that are indicated by the Aw-differences between the 
two sets of data. The cross sectional location of each pair of 
curves is shown by the cross in the inset. Using the outer 
flow solution as the standard for comparison, it is observed 
that for the Aw-differences to vanish, the magnitudes of the 
solid-wall vertical velocities must increase at all field level 
control points except for those that are outboard of the wing 
tip (z/c = 0). Here reductions in the magnitudes of w are 
indicated. The problem is to determine a set of boundary 
conditions on the floor and ceiling panels outboard of the 
wing tip which will increase the magnitudes of the vertical 
velocities at 2 y/b - 1.18 (z/c = 0.52) and which will decrease 
the magnitudes of the velocities at 2 y/b ~ 1 .30 (z/c = 0). This 
problem, which may not exist in an actual wind-tunnel 
experiment, may be associated with the method (Douglas- 
Neumann Program) for computing the flow-field characteris- 
tics near the wing tip. Since the appropriate course of action 
is not known beforehand, the vertical velocities at z/c = 0 
were neglected (for the time being) and the primary consider- 
ation was given to eliminating wall interference effects at the 
other control points. As previously mentioned, an iterative 
procedure would be followed, and at each iteration the 
Aw-differences would be used as inputs to the influence 
coefficient program to determine the wall boundary 
conditions. 

Figure 9 shows the vertical velocity distributions after the 
first iteration. The wall interference was reduced slightly 
with a resulting reduction in Cj^ from 0.218 to 0.204. 

After four iterations a set of vertical velocity distributions 
were obtained showing little or no wall interference at nearly 
all of the field control points. Figure 10 shows the vertical 
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velocity distributions after the four iterations. The data at 
the z/c = ±0.52-control points show little or no indications 
of wall interference. The velocity distributions at the station 
outboard of the tip ( zjc = 0), however, show that wall inter- 
ference still persists there. The lift coefficient of 0.181 is 
higher than the interference-free value of 0.1 61 . 

The convergence rate for this sequence of calculations is 
shown in figure 11. The MSE (mean square error) at each 
iteration is observed to approach, and finally assume, a value 
less than that for the interference-free condition. This is in 
contrast to the MSE of the data outboard of the tip (in the 
wing plane zjc = 0) which becomes larger with each iteration 
(fig. 11(b)). The reason for this is not presently known and 
will be studied further in later investigations. 

Figure 12 shows how C i was reduced with each iteration. 
Here it is evident that although most of the data (at all points 
except at z/c = 0) showed convergence, the did not match 
that for free air. 

Next, an alternate strategy was employed whereby a new 
starting point was found (the second set of iterations). This 
was the result of using boundary conditions which were 
eight times the blowing and suction velocities for the case 
shown in figure 9 ( C L = 0.204). The magnitudes of the 
resulting vertical velocities (fig. 13) are observed to be gener- 
ally greater than those computed using the outer flow solver. 
The most noteworthy result is that for the Aw values to be 
able to vanish, the magnitudes of all the adaptive wall w’s 
must be reduced. (A scale change was required for this figure 
because of the larger magnitudes of the vertical velocities.) 

Seven iterations were required to reach interference-free 
flow conditions. The results of the seventh iteration are 
shown in figure 14. Note that the lift coefficient is close to 
that of the free-air value of 0.161 . 

The convergence rate (fig. 15) shows clearly how the MSE 
for the data at z/c = ±0.52 and those for z/c = 0 progressed 
from a large starting value to a smaller and smaller value with 
each iteration. Although the MSE for the data at the out- 
board station (z/c = 0) did not assume a value less than that 
for free air, it apparently is close enough judging by the close 
agreement of the C L to that of the free-air value of 0.161 . 

Figure 16 shows the effect of iteration on C L . The data 
denoted by the circle symbols are from the previously 
described first set of iterations. The data denoted by the 
square symbols are the results of the second sequence. Both 
the MSE and Ci values show that the adaptive wall wind 
tunnel has effectively achieved free-air conditions. 

Figure 17 shows a comparison between the adaptive wall 
vertical velocities and those for free air. The good agreement 
shown between the two sets of data verifies the adaptive wall 
results. 

To summarize the previously developed procedure, the 
iteration must start from a point where the boundary condi- 
tions for reducing the Aw-differences is consistent from one 
set of field level control points to another. For instance, at 
the control points outboard of the wing tip in figure 8, the 


boundary conditions were to increase the w’s at z/c = ±0.52 
and, simultaneously to decrease the w’s at z/c = 0. This con- 
flicting requirement is not observed in the data in figure 13 
where all the w’s have to be reduced for the Aw values to 
vanish. 

The procedures for establishing interference-free flows 
which were developed for the 2° angle of attack case were 
applied to the wing at a = 1°. Shown in figure 18 are the 
solid wall vertical velocity distributions and those computed 
using the outer flow solver. Using the outer flow solver veloc- 
ities as the standard for comparison, the magnitudes of the 
solid wall velocities must increase at all of the field level con- 
trol points except those at the outboard station (z/c = 0). To 
find the new starting point, the values for Aw were used as 
inputs to the influence coefficient program to solve the 
boundary conditions. These were multiplied by four and 
imposed as boundary conditions in the Douglas-Neumann 
program. The resulting vertical velocity distributions includ- 
ing those from the outer flow solution are shown in fig- 
ure 19. This constituted the starting point from which to 
begin the iteration. Four iterations were required to achieve 
interference-free flow conditions. The results of the final 
iteration are shown in figure 20, where little or no wall inter- 
ference is indicated (except at 2 y/b = 1 .30). 

The rate of convergence is shown in figure 21. The final 
value of the MSE for the data at z/c = ±0.52 is below that 
for free air. The MSE for the data at the points outboard of 
the wing tip (z/c = 0) is shown to approach but not to 
assume values less than that for free air. The lift coefficient 
for the last iteration is 0.086, while that for free air is 0.081. 

Figure 22 shows the effects of iteration for the a = 1° 
case. The a = 2° iterations are also shown in the figure. Note 
that although the C L value was quickly reached, the MSE 
data indicated that three iterations were required to satisfy 
the convergence criteria. 

Adaptive wall vertical velocity distributions and those for 
free air are compared in figure 23. The differences between 
the two sets of data may account for the slight difference in 
C L- 


CONCLUSIONS 


The following conclusions may be drawn from this study: 

1 . The method of influence coefficients may be used to 
obtain a succession of boundary conditions that can be used 
in an iterative scheme to compute interference-free flow in 
the wind tunnel. 

2. The vertical velocities at the upper and lower field level 
control points may be varied independently of each other by 
changing the ceiling and floor boundary conditions, respec- 
tively. This effectively decouples the upper half from the 
lower half of the wind tunnel. 
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3. It is not necessary to introduce blowing/ suction at all 
of the floor and ceiling panels. 

4. The simulation described in this report can be used for 
trade-off studies in the design of three-dimensional adaptive 
wall wind tunnels. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, California, December 15, 1983 
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APPENDIX 


For interference-free flow the vertical velocities computed 
by the outer flow solver would be expected to be identical to 
those computed by the method of reference 24. However, 
since the method of finite differences is used to solve the 
small disturbance equation in the outer flow solver, the accu- 
racy of the solution will depend mainly on the mesh size. In 
the present investigation the mesh was refined several times 
until the Aw-differences (shown in fig. 24) between the free 
air and the outer flow solver data were sufficiently small. 
Because the Aw-differences at the field level control points at 
z/c = ±0.52 were much smaller than those at the z/c = 0- 
control points, two convergence criteria were established 


based on the mean square error. The MSE was specified as 
0.493X10" 3 for the data at z/c = ±0.52 and 1.099X10' 3 for 
the data at z/c = 0. The MSE for the data at z/c = 0 could 
have been reduced further by more reduction in mesh size. 
Since one of the objectives of the present investigation was 
to find a quick method for assessing wall interference, the 
results were considered satisfactory so that further mesh size 
refinement was not required. The criteria adopted herein is 
the following: the flow is considered interference free if the 
MSE of the data at z/c = ±0.52 is less than 0.493X10" 3 and 
if the MSE of the data at z/c = 0 is approximately equal to 
1.099X10” 3 . 
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Figure 2.— Wind-tunnel paneling scheme. 
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NACA 65A006 SECTION 


a L.E. = 3.58° 

c = 8.66 DIMENSIONS IN cm 

Figure 3.- Wing geometry. 



Figure 4.— Coordinate system. 




CEILING PANEL 




Figure 5 Vertical velocity distribution induced by blowing through a single panel point, tunnel empty , M - 0.6. 


WIND TUNNEL 



Figure 6.— Cross sectional view of wind tunnel. 
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(a) 2y/b = 0.24 (b)2y/b = 0.71 


Figure 9 — Vertical velocities for the adaptive wall wind tunnel, first iteration, C, = 0.204. 







(c) 2y/b= 1.18 


Figure 9. 
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x/c 

(d) 2y/b = 1.30 

Concluded. 





(a) 2Y/b = 0.24 (b)2y/b = 0.71 

Figure 10.- Vertical velocities for the adaptive-wall wind tunnel, fourth iteration, C ^ = 0.181. 
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(c) 2y/b =1,18 (d) 2y/b = 1.30 


Figure 10.— Concluded. 
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Figure 1 1 Convergence rate. 
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Figure 12 .— Effect of iteration on C^, a = 2° . 
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(a) CONTROL POINTS 

z/c = 0.52 AND z/c = --0.52 


Figure 15.— 
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ITERATIONS 
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□ 2 nd 



Figure 16.— Effect of iterations on Q- , a = 2°. 
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(a) 2y/b= 0.21 


(b) 2y/b = 0.71 
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z/c = 0.52 


o FREE AIR 
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(c) 2y/b =1.18 (d) 2y/b = 1.30 

Figure 24.— Vertical velocities for interference-free flow, Cj = 0.161 . 
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